Hamiltonian models of multiphoton processes and four-photon 
squeezed states via nonlinear canonical transformations 



Silvio De Siena, Antonio Di Lisi, and Fabrizio Illuminati 
Dipartimento di Fisica, Universita di Salerno, INFM-Unitci di Salerno, and INFN-Sezione di Napoli, 
Gruppo Collegato di Salerno, via S. Allende, 1-84081 Baronissi (Salerno) Italy 

(February 9, 2008) 



O 
O 

00 
(N 



> 



O 
O 



Oh 
I 



We introduce nonlinear canonical transformations that yield ef- 
fective Hamiltonians of multiphoton down conversion processes, and 
we define the associated non-Gaussian multiphoton squeezed states 
as the coherent states of the multiphoton Hamiltonians. We study 
in detail the four-photon processes and the associated non-Gaussian 
four-photon squeezed states. The realization of squeezing, the be- 
havior of the field statistics, and the structure of the phase space dis- 
tributions show that these states realize a natural four-photon gener- 
alization of the two-photon squeezed states. 

PACS numbers: 42.50.Dv, 03.65.-w 

Squeezed states represent a remarkable improvement in 
interferometry, providing an example of clearly nonclassical 
states whose experimental realization is more readily acces- 
sible compared to that, e.g., of the number states. As is well 
known, squeezed states are minimum uncertainty states with 
unbalanced quantum noise on conjugate observables, tuned 
by the squeezing parameter They are typically generated by 
two-photon parametric down conversion processes in media 
with second-order optical nonlinearity. In these systems, at 
variance with the case of one-photon processes, the modes 
of the electromagnetic field are not independently excited but 
instead coupled in pairs. As a consequence, the two-photon 
process is more adequately described by the field quadratures 
rather than by the single field mode operator (the annihilation 
operator) [||]. In recent years, different types of nonclassical 
states and generalizations of the squeezed states have been 
studied intensively in quantum optics as well as in atomic 
physics and in many other physical systems [|[). A chal- 
lenging, yet so far elusive goal has been to find a multipho- 
ton generalization of the effective two-photon Hamiltonian 
model for degenerate down conversion processes Mul- 
tiphoton processes are in fact becoming of current interest in 
the study of the foundations of quantum mechanics. In par- 
ticular, four-photon processes can be implemented by entan- 
gling pairs of two-photon squeezed states in a beam splitter 
[^], while high fidelity teleportation and the realization of EPR 
and GHZ states via four-photon processes have been recently 
demonstrated experimentally [||| . 

In this letter we introduce an effective Hamiltonian model 
of (degenerate) multiphoton down conversion processes and 
multiphoton squeezed states. The model is obtained by gen- 
eralizing the canonical transformations originally exploited in 
[|l|] to define the two-photon squeezed states. The generalized 
transformations achieve the threefold goal of preserving the 
canonical commutation relations, of providing multiphoton 
down-conversion Hamiltonians, and of realizing squeezing. 



We show that these transformations exist and are obtained by 
adding to the linear Bogoliubov transformation arbitrary but 
sufficiently regular nonlinear functions of one of the two field 
quadratures, i.e. the fundamental operators associated to two- 
photon processes 1^. Generalized canonical transformations 
yielding multiphoton processes have been first introduced in 
[0], with the nonlinearity placed on the first field quadrature 
Xi, and the lowest order relevant case of four-photon pro- 
cesses has been discussed as well. However, the photon num- 
ber distribution and the phase space distributions of the asso- 
ciated four-photon squeezed states differ radically from those 
of the two-photon squeezed states already for very small val- 
ues of the strength of the nonlinearity (i.e. the strength of the 
multiphoton contributions). 

Here we define a different class of four-photon squeezed 
states (FPSSs). They are obtained by introducing a nonlinear 
canonical transformation with a quadratic nonlinearity placed 
on the second field quadrature X2, and by antisqueezing it. 
We show that these states, although obtained by a method 
similar to that exploited in [^], are physically very distinct 
from those with nonlinearity on Xi, and are the proper four- 
photon analog of the two-photon squeezed states concerning 
the relevant physical aspects: down conversion Hamiltonian 
and squeezing, field statistics, and phase space distributions. 

Let us consider a single mode a of the electromagnetic 
field. Canonical transformations involving higher powers of 
a and at can be defined, although not directly in terms of 
the field modes, but rather of the field quadratures. There are 
only two possible one-mode generalized canonical transfor- 
mations, and they read: 



bi = + va) + ^lF{Xi) , i = l,2. 



(1) 



Here, a, a' denote the one-mode fundamental canonical vari- 
ables ([a, flt] = 1), bi denotes the transformed mode, F is 
a sufficiently regular, Hermitian function of one of the field 
quadratures Xi — {a + a^)/\/2, X2 = — i(a — at)/\/2, and 
/i, ly, 7 are complex coefficients, where 7 is the "coupling" pa- 
rameter measuring the strength of the nonlinearity. One trans- 
formation is generated by F{Xi), the other by F{X2)- Ac- 
cordingly, the canonical constraint = 1 on the trans- 
formed modes yields two different sets of conditions on the 
coefficients of the transformations. The first condition is com- 
mon to the two sets; it is the standard Bogoliubov constraint 
l/ip — = 1. The second condition, which remarkably 
does not depend on the form of the operatorial function F, is 
either Re{fj,j* — v*"i) ~ 0, or Im{fi'-f* — = 0, depend- 
ing on whether the nonlinearity is placed on Xi or on X2, 
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respectively. Here i?e( ), Im{-) denote the real and the imag- 
inary parts. The two transformations (M yield two different 
one-mode multiphoton Hamiltonians = b\bi + 1/2 and 
Hp = b\b2 + 1/2. It is instructive to write them in terms of 
the field quadratures. Let us denote by Xi the quadrature ar- 
gument of the nonlinear function F, and by Xj the remaining 
quadrature (i, j = 1,2 , i ^ j). Without loss of generality we 
adopt the parameterization /i = coshr , v ~ sinhre^*"^, and 
we choose = 0. We may then write: 



g2n 



-[X^ + V2^,e'''F{X,)]^, 



(2) 



where ri ~ r, and r2 = — r. The coefficient 71 — Im{j), as 
in this case the canonical constraints imply 7 imaginary, while 
72 = 7, as in this case the canonical constraints imply 7 real. 
Finally, Hp ^ H'p if i ^ 1, and Hp = H^p^ if i = 2. It is 
well known that, for the two-photon squeezed state (TPSS), 
the instance r > implies squeezing in Xi as well as anti- 
squeezing in X2, due to the constraint of minimum Heisen- 
berg uncertainty. In the canonical multiphoton extension, this 
is no longer true. In fact, we see from Eq. that the nonlin- 
ear function of Xi is added to the other quadrature Xj. This 
shows that to the same TPSS obtained either by squeezing one 
quadrature or by antisqueezing the other quadrature, there cor- 
respond two different possible, physically distinct multipho- 
ton squeezed states. Squeezing on Xi is implied by the choice 
r > and nonlinearity F onXi, while antisqueezing on X2 is 
implied still by r > but nonlinearity F moved on X2. The 
two distinct states obtained by these two different nonlinear 
transformations reduce to the same TPSS with squeezing on 
Xi (and thus antisqueezing on X2) as the parameter 7 goes to 
zero. The Hamiltonian Hp can be written down completely in 
terms of the fundamental modes a, once a particular form 
of F{X,) is selected. The choice F{X,) = (X^)" , i = 1, 2, 
provides generic 2n-photon Hamiltonians H2,i- The trans- 
formation generated by = (-^^2)^ leads to a four-photon 
Hamiltonian Hl^ = b\b2 + 1/2 that reads: 



H 



II 



3 1 

. — (cosh 2r + 37^)a^a + sinh^r H — 7^^ H — 

4 2 

+ \ie'{a'^ + a) + i(sinh2r - ■i-i'^){a)^ + a^) 

+ ^^2^t2„2 _ ^2(^t3^ + ^t^3) + 1^2(^t4 + „4) 



(3) 



We see that ijj^ contains explicitly two-, three-, and four- 
photon degenerate parametric down conversion terms. The 
multiphoton squeezed states are the eigenvectors {|/3)7.f} 
of the transformed annihilation operator b (here h stands 
generically for either &i or b2), with eigenvalue /3, and 
are the coherent states with respect to the vacuum state 
|0)^.F of the generic multiphoton Hamiltonian: \(3)^^p = 
D{P)\Q)^^p, where D{(3) = e:>qi{0 - (i*b) denotes 
the Glauber unitary displacement operator. The vectors 



{\[3)'y,p} form an overcomplete set and resolve the identity: 
I/tt J d?(3 |/3)-y.F F.i{P\ = I- If, as for the standard squeezed 
states, we choose (3 = /ia + va* , where a = ai + ia2 is 
the coherent amplitude, and again = 0, the canonical trans- 
formations (|l]) are implemented by the unitary operators Ui ~ 
exp [ie''^^iG{Xi)]D{a)S{r). Here S(r) = exp [r[a^^ - a^)] 
denotes the single-mode squeezing operator for = 0, and 
G{x) — F{y)dy. The states \(3)^,p are not minimum un- 
certainty states, but they approximate the minimum uncer- 
tainty on the characteristic time scale f^^ provided by the 
time dependence r(t) = rt of the squeezing parameter [^]. 
In the representation in which the quadrature argument of the 
function F is diagonal, the eigenvalue equation for the states 
\f3)j,p is a simple first order linear differential equation for the 
wave functions 'i>Jj^{xi) = {xi\P)^p , i = 1,2 . Integration 
yields 



^f{x,)^{7:a,)-'/'exp 



Xi X 



(0) 



2ai 



X exp{i[ciXi + V2e'''%G{xi)]} 



(4) 



xf^ 



V2ai, 



-V2, 



a2, ci 



where a-i 

\/2a2, C2 = — \/2cn- The multiphoton squeezed states 
are non-Gaussian because of the non-quadratic term G{xi) 
appearing in the phase of the wave function. However, the 
probability density is Gaussian in the representa- 

tion in which the quadrature argument of the nonlinear func- 
tion F is diagonal. The probability density displays squeezing 
in this quadrature by the usual factor e^'"' . This fact further 
shows the necessary link between squeezing and nonlinear- 
ity for the multiphoton states. From now on we specialize to 
the two possible four-photon cases, respectively associated to 
F = {XiY and to = (^2)^. We work exphcitly in the rep- 
resentation in which the first quadrature Xi is diagonal. In this 
representation we can compare congruently the two different 
four-photon squeezed states (FPSSs), namely '^\{x) associ- 
ated to F = {XiY, and ^'4^(2;) associated to F = {X2Y: 



(ttc ) exp 



— (^a; - V2ai 



X exp 



1V2 



a2X + e ^Iin{'y) — 
o 



*f (x) = exp (fca;)v4i 



Ix + m 



12/3 



(5) 



(6) 



In Eq. (^ iV is the normalization factor, Ai[y] denotes the 
Airy function that goes to zero as y 00, I = e"^ /2j; m = 
e-2'-/i672 _ k = e-747. The states (|) and (|) are 
both non-Gaussian. However, while |^'4(a;)| is Gaussian, 
|\I>4^(x)p is not. Its behavior is shown in Fig. 1. 
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FIG. 1. The probability density \^i'{x)f of tlie FPSS (|) with 
amplitude /3 = 3\/2, squeezing r — 0.8, and coupling 7 — 0.14. 

We now analyze the statistical properties of the state (^). 
In Fig. 2 we plot the photon number distribution P{n) of the 
TPSS (7 = 0) and of the two FPSSs -^i and *f . 
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FIG. 3. The Wigner quasiprobability distribution W{xi,X2) of 
the FPSS (lei) with /3 = 3^2, r = 0.8, and 7 = 0.14. 




FIG. 2. The photon number distribution P{n) at /3 = 3\/2 and 
r = 0.8 for the FPSS ^ with 7 = 0.1 (dot-dashed line); for the 
same state with 7 — 0.5 (doubledot-dashed line); for the FPSS 
with 7 = 0.05 (dashed line); and for the TPSS (full line). 

We see from Fig. 2 that the form of P{n) for the state ^'4^ 
is similar to that of the TPSS and is very stable as a function 
of the nonlinear coupling 7. The only significant difference is 
that, for larger values of the coupling the peak of P{n) moves 
to the left, due to the growing influence of 7 on the mean num- 
ber of photons. We see instead that P(n) for the state is 
unstable and is strongly deformed already for very small val- 
ues of the nonlinear coupling. The profoundly different be- 
havior of P{n) for the states ^'4 and ^4^ is deeply rooted in 
the structure of their respective quasiprobability distributions 
in phase space. In Fig. 3 we plot the Wigner function for the 
state 'i'l^ , while in Fig. 4 we plot it for the state 'i'^. The 
nonclassical nature of the state ^'4^ is more pronounced. Its 
Wigner function is deformed and attains also negative values. 
However, it is neither rotated nor translated with respect to 
that of the TPSS, at variance with the Wigner function of the 
state ^1*4 which is both translated and rotated with respect to 
that of the TPSS. This difference affects crucially the proper- 
ties of phase space interference. The latter is responsible for 
the oscillations of P{n) for the TPSS 




FIG. 4. The Wigner quasiprobability distribution W{xi,X2) of 
the FPSS ^ with l3 = 3^, r = 0.8, and 7 = 0.14. 

The oscillations are preserved in a wide range of values of 7 
for the state ^'4^, exactly because its Wigner function is nei- 
ther translated nor rotated with respect to that of the TPSS, 
while they quickly disappear for the state ^'4, whose Wigner 
function is both translated and rotated with respect to that of 
the TPSS already at small values of 7. In order to complete the 
analysis of the statistical properties of the states ^'4 and ^'4^ 
we move to study their correlation functions. In Fig. 5 we plot 
the normalized second-order correlation functions g'--^^ (0) of 
the TPSS (7 = 0) and of the states and ^'4^, as a function 
of r. We see that the state ^'4^ follows a behavior similar to 
that of the TPSS, being exactly Poissonian at r = 0, slightly 
more sub-Poissonian for r < 0.9, and super-Poissonian for 
larger values of r. The value of saturation of g'^' (0) at asymp- 
totically large values of r is an increasing function of 7. All 
these characteristics remain stable as 7 is varied. At variance 
with this behavior, the state ^'4 is always super-Poissonian, 
even at r = 0, and for small values of 7. 
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FIG. 5. The second-order correlation (?''^'(0) at f3 — 3\/2 as a 
function of r for the four-photon squeezed state ^ with 7 = 0.1 
(doubledot-dashed line); for the FPSS (||) with 7 = 0.1 (dashed 
line); and for the TPSS (full line). 

It is of particular importance to study the fourth-order cor- 
relation function in order to establish the degree of probability 
of simultaneous four-photon detection (bunching) [^. In Fig. 
6 we plot the normalized fourth-order correlation functions 
5(4) (0) of the TPSS (7 = 0) and of the states *i and *f , 
as a function of r. Also in this case the FPSS ^'4^ follows 
a behavior similar to that of the TPSS. It favors four-photon 
anti-bunching for values of r < 0.8. For larger values of r it 
strongly favors four-photon bunching, and the value of satu- 
ration of .g'-"'''(0) grows with 7. The state instead always 
favors four-photon bunching even at 7 = and for small val- 
ues of 7. 
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FIG. 6. The fourth-order correlation g'*'(0) at f3 = 3\/2 as a 
function of r for the FPSS (^ with 7 = 0.1 (doubledot-dashed line); 
for the FPSS (|) with 7 = 0.1 (dashed line); and for the TPSS (full 
line). Inset: the same graph in the interval [0, 1]. 

In summary, we have introduced an effective Hamiltonian 
model for multiphoton down conversion processes via a non- 
linear generalization of the Bogoliubov transformation. This 
generalization allows to define two different classes of non- 
Gaussian multiphoton squeezed states. We have studied in 
detail the four-photon case, and we have shown that the four- 
photon squeezed states associated to quadratic nonlinearity in 
the second field quadrature and to antisqueezing in the same 



quadrature are the natural four-photon extension of the two- 
photon squeezed states with respect to the Hamiltonian struc- 
ture, the field statistics, and the Wigner quasiprobability dis- 
tribution. We have considered the instance of degenerate mul- 
tiphoton parametric down conversions. Work is in progress 
to address the description of nondegenerate, two- and multi- 
mode multiphoton processes. 
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